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Recent  experiments  have  suggested  that  there  is  a  critical  size  for  patterned  silicon  (Si)  thin  film  electrodes 
for  delamination  from  a  current  collector  during  lithiation  and  delithiation  cycling.  However,  no  existing 
theories  can  explain  this  phenomenon,  in  spite  of  its  potential  importance  in  designing  reliable  electrodes 
for  high-capacity  lithium-ion  batteries.  In  this  study,  we  show  that  the  observed  delamination  size  effect 
can  be  rationalized  by  modeling  thin  film  delamination  in  the  presence  of  large  scale  interfacial  sliding. 
A  method  is  proposed  to  deduce  the  critical  size  for  delamination  based  on  the  critical  conditions  for 
the  nucleation  and  growth  of  edge  or  center  cracks  at  the  film-substrate  interface  under  plane  strain 
or  axisymmetric  conditions.  Applications  to  lithiation  of  thin-film  Si  islands  give  results  in  excellent 
agreement  with  experimental  observations. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Mechanical  degradation  such  as  fracture  and  delamination  of 
electrodes  due  to  lithiation  and  delithiation  has  been  widely  rec¬ 
ognized  as  one  of  the  main  causes  for  the  short  cycle  life  of 
high-capacity  lithium-ion  (Li-ion)  batteries  [1,2].  For  such  appli¬ 
cations,  silicon  (Si)  holds  the  maximum  theoretical  charge  capacity 
and  has  been  intensely  studied  as  one  of  the  most  promising 
negative  electrode  materials.  However,  mechanical  degradation 
has  been  frequently  reported  in  Si  electrodes  fabricated  in  var¬ 
ious  forms  of  thin  films  [2-7],  nanowires  [8,9],  and  composite 
nano-structures  [10].  It  has  been  shown  that  250  nm  thick  Si 
films  on  copper  (Cu)  foils  exhibit  near  theoretical  charge  capac¬ 
ity  for  a  limited  number  of  cycles,  but  the  capacity  fades  drastically 
upon  further  cycling  due  to  fracture  and  delamination  [7].  Numer¬ 
ous  studies  have  been  devoted  to  improving  the  electrochemical 
and  mechanical  stability/performance  of  Si  electrodes  by  reduc¬ 
ing  the  size  of  Si  particles  below  a  critical  fracture-resistant  size 
[6],  enhancing  adhesion  between  thin  film  electrodes  and  substrate 
[11],  avoiding  formation  of  crystalline  Lii5Si4  phase  which  is  associ¬ 
ated  with  formation  of  high  internal  stresses  [12,13],  and  reducing 
the  film  thickness  to  achieve  higher  cycling  rates  [14].  For  thin- 
film  Si  electrodes  deposited  on  substrates,  it  is  worth  pointing  out 
that,  although  the  lithiation  and  delithiation  cycling  may  generate 
microcracks  through  the  film  thickness,  the  fractured  film  can  still 
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serve  as  the  active  material  as  long  as  it  can  remain  adhered  to  the 
current  collector.  In  fact,  through-the-thickness  cracks  in  such  sys¬ 
tems  can  alleviate  diffusion-induced  stresses  by  facilitating  volume 
changes  associated  with  Li  insertion/extraction.  Therefore,  while 
fracture  in  thin  film  electrodes  partly  decreases  their  capacity  due 
to  the  formation  of  solid-electrolyte  interphase  (SEI)  layers  at  newly 
created  cracked  surfaces,  film  delamination  is  a  far  more  critical 
issue  to  capacity  fading  in  such  electrodes. 

The  present  paper  is  aimed  to  present  a  method  to  deduce  the 
critical  size  for  interfacial  delamination  of  thin  film  electrodes  by 
examining  the  critical  conditions  for  the  nucleation  and  propa¬ 
gation  of  interfacial  cracks  at  the  film-substrate  interface  under 
plane  strain  or  axisymmetric  conditions.  The  motivation  for  this 
study  has  come  from  the  following  experimental  observations  on 
the  delamination  behaviors  of  Si  thin  film  electrodes  which  have 
been  reported  in  our  recent  paper  [15]  but  remains  unexplained  to 
this  date.  First,  it  was  observed  that  a  continuous  Si  thin  film  with 
thickness  of  lOOnm  exhibits  through-the-thickness  cracks  with 
spacing  in  the  range  of  5-10  pim  due  to  lithiation  and  delithiation. 
Next,  a  set  of  experiments  was  performed  to  investigate  the  frac¬ 
ture  behaviors  of  patterned  Si  thin  film  square  islands  on  Cu  with 
the  same  film  thickness  of  lOOnm  but  3  different  lateral  dimen¬ 
sions:  7  pun  x  7  pim,  17  pim  x  17  pim  and  40  pim  x  40  pim.  Scanning 
electron  microscope  images  then  showed  that  the  7  pim  x  7  pim  Si 
islands  stayed  well  adhered  to  the  substrate,  while  a  fraction  of 
the  17  pim  x  17  pim  islands  and  all  of  the  40  pim  x  40  pim  islands 
peeled  off  the  substrate  after  30  charge-discharge  cycles.  There¬ 
fore,  a  critical  island  size  seems  to  exist  in  the  order  of  10-20  pim 
for  interfacial  delamination  of  patterned  thin  films  with  thickness 


358 


H.  Haftbaradaran  et  al.  /  Journal  of  Power  Sources  206  (2012)  357-366 


bU  vU  um 

□  (TOC 


Fig.  1.  Micrographs  demonstrating  lateral  expansion  and  contraction  of  Si  islands 
on  Ti  current  collector  due  to  lithiation  and  delithaition.  (a)  As-deposited 
17  |jim  x  17  |jim  Si  islands;  (b)  lateral  dimensions  of  the  Si  islands  expand  by  about 
30%  parallel  to  the  substrate,  demonstrating  sliding  along  the  interface;  (c)  Si  islands 
return  to  their  original  size  after  delithiation.  The  scale  bar  is  50  pirn  (after  [21  ]). 


of  1 00  nm.  The  existing  theories  of  interfacial  delamination  [  1 6-20] 
simply  cannot  explain  why  the  lOOnm  thick  islands  exhibit  such 
size-dependent  delamination  behavior  in  the  size  range  nearly  two 
orders  of  magnitude  beyond  the  film  thickness.  In  the  present  study, 
we  show  that  this  effect  can  only  be  rationalized  by  considering  thin 
film  delamination  in  the  presence  of  large  scale  interfacial  sliding. 
Fig.  la  shows  an  array  of  17  p>m  x  17  p>m  Si  islands  deposited  on 
a  titanium  (Ti)  current  collector.  Upon  lithiation,  it  was  observed 
that  the  Si  islands  underwent  about  30%  of  lateral  expansion  par¬ 
allel  to  the  substrate,  as  demonstrated  in  Fig.  lb.  Fig.  lc  shows  that 
the  same  islands  were  able  to  shrink  back  to  their  original  size  upon 


deilthiation.  These  in-plane  expansion  and  contraction  of  Si  islands 
indicate  large  scale  sliding  along  the  film-substrate  interface  [21]. 
By  accounting  for  the  effect  of  interfacial  sliding,  we  will  develop 
a  method  to  deduce  the  critical  size  for  interfacial  delamination 
in  patterned  electrode  structures  by  examining  the  critical  condi¬ 
tions  for  the  nucleation  and  propagation  of  edge  or  center  cracks 
at  the  film-substrate  interface  under  plane  strain  or  axisymmetric 
conditions. 

It  is  worth  pointing  out  that  the  present  work  can  be  regarded 
as  part  of  a  series  of  studies  on  patterned  thin  film  electrodes 
[15,21,22].  In  a  previous  study  [15],  we  have  already  proposed 
a  shear  lag  model  which,  by  accounting  for  interfacial  sliding  or 
substrate  yielding,  explains  the  experimentally  observed  fracture 
spacing  in  continuous  Si  thin  films.  This  model,  when  applied  to 
patterned  Si  thin  film  islands,  would  imply  that  there  is  substan¬ 
tial  sliding  at  the  film-substrate  interface  and  that  it  is  possible 
to  choose  the  island  size  such  that  the  maximum  stress  in  the 
island  stays  below  the  flow  stress  of  lithiated  Si,  in  which  case 
the  island  would  remain  elastic  during  Li  insertion/extraction  [22]. 
The  same  model  has  been  successfully  used  to  model  stress  evo¬ 
lution  in  patterned  thin  film  islands  on  Ti  substrates,  with  results 
in  excellent  agreement  with  experiments  [22].  According  to  these 
experiments,  the  interfacial  sliding  strength  has  been  estimated  at 
about  40  MPa  between  lithiated  Si  and  Cu  [15],  and  about  10  MPa 
between  lithiated  Si  and  Ti  [21,22],  which  are  comparable  with 
direct  experimental  measurements  of  interfacial  sliding  resistance 
in  other  similar  kinds  of  material  systems  [23].  In  spite  of  these 
progresses,  the  anomalous  critical  size  for  interfacial  delamination 
in  thin-film  Si  electrodes  during  lithiation  and  delithiation  has  not 
been  explained  and,  therefore,  the  question  of  what  factors  con¬ 
trol  the  delamination  of  electrodes  from  the  current  collector  still 
remains  open.  In  the  present  manuscript,  we  shall  attempt  to  close 
this  gap  by  investigating  the  delamination  behaviors  of  patterned 
thin  film  electrodes  in  the  presence  of  large  scale  interfacial  sliding. 


2.  Theories  of  fracture  and  delamination 

Theories  of  fracture  mechanics  and  diffusion-induced  stresses 
have  been  employed  to  model  fracture  of  electrodes,  with  results 
suggesting  the  existence  of  a  critical  size  and  charging  rate  below 
which  fracture  becomes  impossible  [6,24-31].  A  widely  accepted 
result  is  that  sufficiently  thin  electrodes  are  resistant  against  frac¬ 
ture  in  the  thickness  direction. 

Delamination  of  thin  films  on  substrate  has  been  widely  stud¬ 
ied  in  the  literature  [16-20].  An  assumption  commonly  made  in 
modeling  delamination  of  thin  films  on  substrate  is  that  the  film  is 
well-bonded  to  the  substrates  with  no  interfacial  sliding  [16-20]. 
Under  plane  strain  conditions,  a  well-known  result  is  that  the 
energy  release  rate  associated  with  an  edge  delamination  crack 
reaches  a  steady-state  value  of  G0  =  (1  -  v2)cr2h/2E  as  soon  as  the 
crack  size  becomes  comparable  to  or  larger  than  the  film  thickness 
ft,  where  ay  is  the  stress  in  the  film,  and  E  and  v  are  the  Young’s  mod¬ 
ulus  and  Poisson’s  ratio  of  the  film,  respectively  [16,19].  However, 
when  the  interfacial  sliding  strength  is  much  smaller  than  the  stress 
in  the  film,  as  is  the  case  expected  for  a  Si  electrode  on  a  current 
collector,  a  large  scale  interfacial  sliding  zone  would  be  expected 
near  the  edges  of  the  film.  If  the  length  of  the  film  is  larger  than 
the  size  of  the  sliding  zone,  a  delamination  crack  would  grow  in  a 
steady  state  manner  as  soon  as  the  crack  size  is  larger  than  the  film 
thickness.  In  this  state,  the  energy  release  rate  is  independent  of  the 
crack  size,  and  there  is  no  size  effect  (Fig.  2a).  On  the  other  hand,  if 
the  interfacial  sliding  zones  from  the  island  edges  are  so  large  that 
they  extend  over  the  entire  island  (Fig.  2b),  a  strong  delamination 
size  effect  would  emerge  over  the  size  scale  associated  with  the  size 
of  the  sliding  zones,  rather  than  the  film  thickness.  The  shaded  area 
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Fig.  2.  Schematic  of  an  edge  delamination  crack  in  a  thin  film  on  substrate  in  the 
presence  of  large  scale  interfacial  sliding,  (a)  The  crack  grows  in  a  steady  state  man¬ 
ner  when  the  film  width  is  much  larger  than  the  interfacial  sliding  zone  at  the  film 
edge,  (b)  The  interfacial  delamination  becomes  dependent  on  the  island  size  as  the 
sliding  zones  from  island  edges  overlap.  In  this  case,  the  energy  release  rate  for 
interfacial  delamination  is  expected  to  decrease  as  the  island  size  is  reduced. 


in  Fig.  2b  shows  the  region  in  which  stress  is  released  as  two  edge 
cracks  nucleate  and  grow  to  a  size  a.  It  is  evident  that  the  shaded 
area  in  Fig.  2b  would  shrink  as  the  island  size  is  reduced,  and  hence 
the  energy  release  rate  is  expected  to  decrease.  This  is  in  contrast 
to  the  case  shown  in  Fig.  2a,  which  corresponds  to  a  relatively  large 
island  size  compared  to  the  size  of  the  sliding  zones,  leading  to  a 
size-independent  energy  release  rate.  Since  the  size  of  the  inter¬ 
facial  sliding  zones  also  determines  the  length  scale  over  which 
the  normal  stress  in  the  film  reaches  its  maximum  value,  the  size 
effects  associated  with  interfacial  sliding  immediately  suggests  an 
explanation  for  the  experimentally  observed  correlation  between 
the  fracture  spacing  in  a  continuous  Si  thin  film  and  the  critical  size 
for  interfacial  delamination  in  patterned  Si  islands.  Therefore,  inter¬ 
facial  sliding  between  lithiated  Si  electrodes  and  substrate,  which 
can  occur  on  a  much  larger  length  scale  than  the  film  thickness, 
is  a  key  to  deducing  the  critical  size  for  interfacial  delamination 
in  patterned  Si  electrodes.  The  reason  that  the  existing  theories  of 
interfacial  delamination  [16-20]  cannot  explain  the  experimental 
observation  that  Si  islands  with  thickness  of  lOOnm  and  length 
of  7  p>m  stayed  well  adhered  on  the  substrate  while  those  with 
the  same  thickness  but  length  of  40  p>m  delaminated  away  can 
be  attributed  to  the  fact  that  the  assumption  of  well-bonded  and 
non-slipping  interface,  while  applicable  to  many  conventional  thin 
film-substrate  systems,  no  longer  holds  in  the  case  of  lithiation  of 
high-capacity  electrodes. 

In  the  following,  we  will  deduce  the  critical  size  for  interfacial 
delamination  by  detailed  plane  strain  and  axisymmetric  analysis 
of  potential  nucleation  and  growth  of  edge  and  center  interfacial 
cracks  in  the  presence  of  large  scale  interfacial  sliding.  We  will  take 
advantage  of  the  symmetry  in  the  problem  to  simplify  the  analysis. 
Furthermore,  we  will  treat  the  edge  and  center  cracks  separately, 
since  their  coexistence  would  decrease  the  elastic  energy  release 
rate  in  the  islands  as  the  driving  force  for  crack  nucleation  and 


Fig.  3.  Schematic  of  (a)  a  Si  island  on  substrate  with  edge  delamination  cracks  and 
(b)  a  free-body  diagram  of  the  island. 


growth.  The  interface  between  the  island  and  substrate  is  assumed 
to  be  well-bonded  under  small  shear  tractions  but  interfacial  slid¬ 
ing  is  allowed  as  soon  as  the  interfacial  shear  stress  reaches  a  critical 
value  r0.  It  is  also  assumed  that  the  sliding  follows  a  perfectly  plas¬ 
tic  response  so  that  the  interfacial  shear  stress  never  exceeds  the 
sliding  strength  r0.  The  immediate  consequence  of  this  behavior  is 
that  the  interfacial  sliding  zones  start  to  develop  at  the  edges  of  the 
film  and  spread  toward  the  center  of  the  island.  Once  they  spread 
to  the  center  of  the  island,  the  entire  island  would  be  sitting  on  top 
of  a  slippery  surface  and  can  freely  expand  upon  further  lithiation. 
This  is  the  case  as  long  as  deformation  within  the  island  remains 
elastic,  i.e.  the  maximum  stress  in  the  island  does  not  exceed  the 
flow  stress  of  lithiated  Si. 


3.  Nucleation  and  growth  of  edge  delamination  cracks  in 
the  presence  of  large  scale  interfacial  sliding 

In  this  section,  we  will  investigate  the  nucleation  and  growth  of 
edge  delamination  cracks  in  a  Si  island  with  interfacial  sliding  zones 
spreading  over  the  entire  island.  For  simplicity,  it  will  be  assumed 
that  the  film  is  thin  enough  that  variation  of  stress  in  the  thickness 
direction  can  be  neglected. 


3. 1 .  Plane  strain  analysis 


Consider  a  Si  island  of  length  21  and  thickness  h  on  a  substrate, 
with  edge  delamination  cracks  of  size  a  at  the  interface  at  both 
edges  of  the  island,  as  shown  in  Fig.  3.  Assuming  that  the  interfacial 
sliding  zones  from  both  edges  of  the  island  overlap,  the  axial  stress 
distribution  in  half  of  the  patch  would  be  given  as 

_r0J/-a-x/-a>x>0  m 

axx-T\0  l  >  x  >  l  -  a  ’  UJ 


where  x  is  distance  from  the  center  of  island.  The  total  elastic  energy 
stored,  per  unit  out-of-plane  length,  in  half  of  the  island  is  obtained 
as 


>+/ 


Ud  = 


1  -  V' 

2  E 


a 


n  (1  -  v2)r?  1  o 

2  dx=  — — -  -(/  -  a)3, 


XX 


2hE 


(2) 


where  E  and  v  are,  respectively,  the  Young’s  modulus  and  Poisson’s 
ratio  of  the  film.  To  investigate  the  criterion  for  crack  growth,  we 
first  note  that  the  elastic  energy  release  rate  G  is 


dUei 

da 


(1  -  V2)l -2 
2hE 


(l -a)2. 


(3) 


Interface  delamination  is  accompanied  with  formation  of  new  sur¬ 
faces  which  increases  the  total  free  energy  of  the  system.  Denoting 
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the  interface  fracture  energy  (toughness)  by  r,  the  total  surface 
energy  of  the  system  Usurf  increases  at  the  rate  of 


dUsurf 
~dcT  =  F 


(4) 


The  Griffith  condition  for  crack  growth  is  that  the  elastic  energy 
release  rate  must  exceed  the  growth  rate  of  surface  energy,  that  is 
G  >  dUsurflda ,  which  yields 


( l-a)> 


2hEf 


(1  -  v2)z 2 


(5) 


For  an  initially  crack-free  interface,  i.e.  a  =  0,  Eq.  (5)  suggests  that 
edge  delamination  at  the  interface  is  only  favorable  for  islands 
larger  than  a  critical  length, 


/  >  / 


2hEf 


cr,nuc  — 


(1  -  v2)t 2 


(6) 


In  this  case,  nucleation  of  edge  interfacial  cracks  decreases 
the  total  free  energy  of  the  system.  Fig.  4a  shows  the  varia¬ 
tions  of  the  elastic  energy  rate  vs.  the  crack  length  a,  where  the 
dimensionless  parameters  Uei  =  2/iFHe//(l  -  v2)r2l3  and  Usurf  = 
2hEUsurf/(l  -  v2)Xq13  have  been  introduced.  As  shown  in  this  fig¬ 
ure,  the  elastic  energy  release  rate  associated  with  crack  growth 
decreases  with  the  crack  length,  while  surface  energy  rises  at  a 
constant  rate.  Therefore,  the  crack  is  expected  to  be  arrested  at  a 
residual  length  aQrr.  This  residual  length  can  be  obtained  from  the 
condition  that  the  total  released  elastic  energy  of  the  system  prior 
to  arresting  has  been  completely  consumed  to  create  new  surfaces, 
and  thus  there  is  no  more  energy  available  for  crack  to  grow,  i.e. 


Uel(a  =  0)  -  Uel(a  =  dan)  =  Usurf(a  =  aan)  -  Usurf(a  =  0).  (7) 

Straightforward  manipulation  shows  that 

aarr  =  2^  ~  \J cr,nuc  ~  f2))-  (8) 

Fig.  4b  shows  how  aarr//  varies  with  l\lcr,nuc •  One  can  easily  verify 
that  at  this  crack  length,  the  elastic  energy  release  rate  becomes  less 
than  the  surface  energy  growth  rate.  An  interesting  case  happens 
when  ciarr  >  l  corresponding  to  the  scenario  when  the  two  edge 
cracks  would  grow  all  the  way  to  the  center  of  the  island.  In  this 
case,  the  island  will  peel  off  from  the  substrate.  It  can  be  shown  that 
this  is  the  case  only  when  the  length  of  the  island  is  larger  than  a 
critical  length  associated  with  full  delamination, 


^  —  l cr,del 


=  73/ 


6hEf 


cr,nuc  — 


(1-v2)t2 


(9) 


To  further  elucidate  the  delamination  mechanisms,  Fig.  4c  shows 
changes  in  the  dimensionless  elastic  and  surface  energies  of  the  sys¬ 
tem,  i.e.  | AGe/|  =  Uei( 0)  —  Uei(a ),  and  ALfsur^-  =  Hsurj(cz)  —  USUrf( 0), 
for  various  island  sizes,  as  a  function  of  the  crack  length.  It  can 
be  seen  that,  in  the  critical  case  of  l  =  lcr  dei  =  V3lcr,nuc,  marked  by  a 
solid  circle  in  Fig.  4c,  the  total  released  elastic  energy  is  just  enough 
to  drive  the  crack  all  the  way  to  the  center.  As  is  evident  from  the 
same  figure,  larger  islands  with  l  =  2lcr,nuc>kr,dei  have  even  higher 
elastic  energy  and  will  delaminate  away.  However,  for  a  shorter 
island  with  l  =  1 .5 lcr,nuc <  lCr,deb  the  surface  energy  grows  beyond  the 
total  elastic  energy  available  in  the  island,  indicating  that  the  full 
delamination  is  impossible,  in  consistence  with  Eq.  (9).  In  the  latter 
case,  the  ultimate  length  of  the  arrested  crack  can  be  obtained  via 
Eq.  (8). 

A  couple  of  remarks  are  worth  making  at  this  point.  While  in 
the  above  analysis  we  have  considered  a  pair  of  symmetrical  edge 
interfacial  cracks,  one  can  verify  through  straightforward  analy¬ 
sis  that  consideration  of  a  single  edge  crack  does  not  change  the 
results.  Furthermore,  in  the  present  analysis,  we  have  ignored  the 


cr.nuc 


Fig.  4.  Plane  strain  analysis  of  interfacial  delamination  of  a  Si  island  on  substrate, 
(a)  Elastic  energy  release  rate  and  surface  energy  growth  rate  vs.  crack  length  for 
various  island  sizes;  (b)  length  of  the  arrested  edge  crack  as  a  function  of  the  island 
size;  (c)  changes  in  the  elastic  energy  and  surface  energy  as  a  function  of  the  crack 
length  for  different  island  sizes. 


elastic  energy  stored  in  the  substrate.  Treating  substrate  as  an  elas¬ 
tic  half-plane,  a  dimensional  analysis  would  reveal  that  the  ratio 
of  elastic  energy  stored  in  the  substrate  to  that  in  the  film  scales 
as  F/i/Fs(/-a),  where  Es  is  the  Young’s  modulus  of  the  substrate. 
Therefore,  while  our  analysis  is  exact  for  the  case  of  a  rigid  sub¬ 
strate,  it  remains  valid  even  for  an  elastic  substrate  as  long  as  the 
length  of  the  island  in  contact  with  the  substrate  is  much  larger  than 
the  thickness  of  the  island.  This  is  indeed  the  case  for  patterned  Si 
thin  film  electrodes  with  length  on  the  order  of  tens  of  microm¬ 
eters  and  thickness  of  a  few  hundred  nanometers.  As  the  island 
loses  its  contact  with  the  substrate  along  the  interface  with  con¬ 
tact  area  shrinking  to  a  size  comparable  to  the  film  thickness,  the 
above  assumption  starts  to  break  down.  However,  when  this  hap¬ 
pens,  the  island  can  be  considered  as  having  effectively  detached 
from  the  substrate. 

For  further  physical  insights,  some  numerical  estimates  of  the 
critical  length  scales  for  interfacial  delamination  would  be  helpful. 
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CT«,  =  I¥((2  +  v)(R_a)_(1+2v)r)'  (11) 

For  the  stress  field  given  in  Eqs.  (10)  and  (11),  the  strain  energy 
density  u,  within  the  region  R  -  a  >  r  >  0,  is  calculated  as1 

U  =  T°(81£~2V)(2 *(2  +  V)2(R  -  a)2  -  6(1  +  v)(2  +  v)r(R  -  a) 

+  (1  +  v)(5+4v)r2).  (12) 


The  total  elastic  energy  stored  in  the  film  is  obtained  by  integrating 
u  over  the  entire  island,  i.e.  Ud  =  fQ  2nrhu  dr,  which  yields 


Ua  =  *Jl^±2>^-af. 


(13) 


The  elastic  energy  release  rate  is  then 


dUel  ;r(l  -  v)(5  +  v)  rl 
=  ~ib= - 9 - m{R~a) 


(14) 


The  formation  of  the  annular  edge  crack  causes  an  increase  in  sur¬ 
face  energy, 

Usurf  =  7T(R2-(R-a)2),  (15) 


which  grows  at  a  rate  of 

=  2jt{R  -  a)r .  (16) 

Using  Eqs.  (14)  and  ( 1 6),  the  Griffith  condition  for  crack  growth,  i.e. 
G  >  d\Jsurfjda,  gives 


(R  -  a)  > 


18  Ehr 


(1  -  v){5  +  v)r2 


2  ' 


(17) 


shear  stress 

Fig.  5.  Schematic  of  (a)  a  Si  circular  island  on  substrate  with  an  annular  edge  crack 
and  (b)  a  free-body  diagram  of  the  island. 


For  a  lOOnm  thick  Si  patch  with  £  =  40  GPa  and  v  =  0.22  (for  lithi- 
ated  Si  [32]),  r0  ~  10  MPa  (for  patterned  Si  thin  films  [21,22]),  and  a 
typical  value  of  1  J/m2  for  the  interface  energy  r,  the  critical  island 
sizes  for  crack  nucleation  and  full  delamination  are  estimated  to 
be  2lCr,nuc^  18.3  |Jtm  and  2/crde/^31  p>m,  respectively,  correspond¬ 
ing  to  maximum  stresses  of  0.92  GPa  and  1.59  GPa  in  the  center 
of  the  island  prior  to  delamination.  These  results  are  in  excellent 
agreement  with  the  experimentally  observed  critical  island  sizes 
for  delamination  and  the  typical  values  of  film  stress  when  this 
happens. 

3.2.  Axisymmetric  analysis 

In  this  section,  we  consider  a  circular  island  on  a  flat  substrate 
and  investigate  the  criterion  for  the  nucleation  and  growth  of  an 
edge  delamination  crack  under  axisymmetric  conditions.  Following 
a  similar  approach  adopted  in  the  preceding  section,  we  assume 
that  the  intercalation  strain  is  large  enough  to  cause  sliding  over  the 
entire  island.  The  radius  of  the  island  is  denoted  by  R  and  the  edge 
delamination  crack  considered  here  is  an  annular  concentric  ring 
whose  front  is  located  at  distance  a  from  the  edge,  as  shown  in  Fig.  5. 
Treating  the  axisymmetric  interfacial  shear  traction  as  a  constant 
radial  body  force  in  the  island  and  assuming  that  the  deformation  of 
the  island  can  be  described  within  plane  stress  approximation,  the 
radial  and  tangential  components  of  stresses  in  the  film,  at  distance 
r  from  the  center  of  island,  can  be  obtained  as  (see  Appendix  A). 

grr  =  ^3  V^y(R-Q-r),  (10) 


suggesting  that  crack  nucleation  in  an  initially  crack-free  interface 
is  favorable  only  for  islands  larger  than  a  critical  size, 


R  >  R 


18  Ehr 


cr,nuc  — 


(1- v)(5  +  v)rg 


2  * 


Similar  to  Eq.  (7),  the  crack  would  arrest  when 


(18) 


aarr  —  R  ~  yj 2R cr,nuc  ~  ^2-  (19) 

Note  this  equation  is  valid  as  long  as  R  <  V2Rcr,nuc •  For 

R  ^  Rcr,del  =  V^R cr,nuct  (20) 

crack  arrest  is  impossible  and  the  edge  crack  would  grow  all  the 
way  to  the  center,  leading  to  complete  delamination. 

For  a  1 00  nm  thick  Si  circular  patch  with  E  =  40  GPa  and  v  =  0.22 
[32],  ro~10MPa  [21,22]  and  r  =  lj/m2,  numerical  estimates  for 
the  critical  island  sizes  for  nucleation  and  delamination  by  axisym¬ 
metric  delamination  are  2Rcrmc  ~  26.6  p>m  and  2Rcr  dei  ^37.6  |jim, 
corresponding  to  maximum  radial  stresses  of  0.98  and  1.39  GPa  at 
the  center  of  the  island.  Fig.  6a  shows  variations  of  the  dimension¬ 
less  elastic  energy  release  rate  and  surface  energy  growth  rate  as 
a  function  of  crack  size.  In  this  figure,  Ud  =  18/iELle//7r(l  -  v)(5  + 

v)TqR4  and  Usurf  =  18/iELlsur^/7r(l  -  v)(5  +  v)TqR4  are  dimension¬ 
less  elastic  and  surface  energies,  respectively.  As  shown  in  this 
figure,  at  a  =  0,  corresponding  to  the  case  of  no  preexisting  crack, 
the  elastic  energy  release  rate  is  equal  to  the  rate  of  increase  in 


1  In  the  analysis  presented  here,  we  have  neglected  the  elastic  energy  stored  in 

the  delaminated  part  of  the  film  for  two  reasons.  First,  we  have  carried  out  the  full 

analysis  that  includes  this  part  of  the  elastic  energy  and  found  no  substantial  changes 

in  the  results.  Second,  a  delaminated  thin  membrane  is  susceptible  to  wrinkling 
which  would  cause  part  of  the  elastic  energy  to  be  released,  bringing  the  solution 
closer  to  the  present  simplified  analysis. 
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Fig.  6.  Axisymmetric  analysis  of  interfacial  delamination  of  a  circular  Si  island  on 
substrate,  (a)  Elastic  energy  release  rate  and  surface  energy  growth  rate  vs.  crack 
length  for  various  island  sizes;  (b)  length  of  the  arrested  edge  crack  as  a  function  of 
the  island  size;  (c)  changes  in  the  elastic  energy  and  surface  energy  as  a  function  of 
the  crack  length  for  different  island  sizes. 


surface  energy  when  R  =  Rcr,nuc •  When  R  >  Rcr,nuc  {R  <  Rcr,nuc ),  the  elas¬ 
tic  energy  release  rate  becomes  larger  (smaller)  than  the  surface 
energy  growth  rate  at  a  =  0,  indicating  that  crack  nucleation  is  ener¬ 
getically  favorable  (unfavorable),  since  it  decreases  (increases)  the 
total  free  energy  of  the  system. 

Fig.  6b  plots  the  dimensionless  residual  length  of  the  arrested 
crack,  aarr/ft,  as  a  function  of  the  dimensionless  radius  of  the  island 
R/Rcr,nuc •  It  can  be  seen  from  this  figure  that  the  length  of  the 
arrested  crack  would  reach  the  radius  of  the  island  (i.e.aarr  =  R) 
when  R  =  Rcr,dei  =  V2 Rcr,nuc ,  corresponding  to  full  delamination. 
Fig.  6c  plots  the  changes  in  dimensionless  elastic  and  surface  ener¬ 
gies,  |  AL/ej|  =  Uei( 0)  —  Uei(a)  and  A 0SUj-f  =  Usurf(a)  ~  ^surf(^X  for 
islands  of  different  sizes  as  a  function  of  the  crack  length.  It  indicates 
that  for  an  island  with  R  =  Rcr,dei  =  y/2RCr,nuc,  the  total  released  elas¬ 
tic  energy  is  just  sufficient  to  drive  the  crack  all  the  way  to  the 
center.  This  critical  case  has  been  marked  by  a  solid  circle  in  Fig.  6c. 
Further  discussions  are  similar  to  the  plane  strain  analysis. 


2a 


Fig.  7.  Schematic  of  (a)  a  Si  island  on  substrate  with  a  center  crack  and  (b)  a  free-body 
diagram  of  the  island. 


4.  Analysis  of  center  delamination  cracks  in  the  presence  of 
large  scale  interfacial  sliding 

In  this  section,  we  consider  interfacial  cracks  at  the  center  of 
islands  using  a  similar  energy  approach.  Under  the  same  assump¬ 
tions  as  stated  in  Section  3,  we  consider  the  analysis  under  plane 
strain  and  axisymmetric  conditions  separately. 


4.2.  Plane  strain  analysis 


Consider  the  same  plane  strain  problem  as  in  Section  3.1,  but 
with  an  interface  crack  of  length  2 a  located  at  the  center  of  the 
island,  as  shown  in  Fig.  7.  Under  the  same  assumption  of  over¬ 
lapping  interfacial  sliding  zones,  the  normal  stress  in  the  film  at 
distance  x  from  the  center  of  island  can  be  written  as 


<3 xx 


To 

h 


l  -  a  a  >  x  >  0 
l-x  l >  x > a 


(21) 


It  follows  that  the  total  elastic  energy  in  half  of  the  island  is 


d 


Uei  = 


(1  _  V2)  n  (1  -  vl)xi  n 

^  olv  dx  = - — — -(/  -  a)2(/  +  2a), 


o 


2  E  6 Eh 

and  the  elastic  energy  release  rate  as 
dUel  (l-v2)r2 


(22) 


G  = 


da 


HE 


a(l  -  a). 


(23) 


Note  that  the  energy  release  rate  vanishes  at  a  =  0,  indicating  that 
a  preexisting  crack  is  needed  for  the  crack  growth  condition  G>  T. 
This  is  in  sharp  contrast  to  the  case  of  edge  delamination  where  no 
preexisting  cracks  are  necessary.  As  shown  in  the  preceding  section, 
even  in  an  initially  crack-free  island,  it  is  possible  for  an  edge  delam¬ 
ination  crack  to  nucleate  spontaneously.  However,  this  is  not  the 
case  for  center  cracks.  To  satisfy  the  growth  condition,  the  length 
of  the  preexisting  center  crack  a0  should  fall  in  the  following  range 


Qo  ,max  >  ao  >  ^0,min* 

where 


ao. 


min 


=  1 


1 


4  FEh 


^0,max  —  2  ^  1  T 


1  - 


(1  -  v2)rp2 


4  FEh 


(1  -  v2)r2/2 

The  above  equation  breaks  down  for 


r> 


(1  -  V2)/2Tq 


and 


4  Eh 


or  l  <  l*r  = 


AEhf 


(1  -v2)r2' 


(24) 


(25) 
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Fig.  8.  Plane  strain  analysis  of  interfacial  delamination  of  a  Si  island  on  substrate  with  a  center  crack  at  the  interface:  (a)  elastic  energy  release  rate  and  surface  energy  growth 
rate  vs.  crack  length  for  various  island  sizes;  total  elastic  energy  of  the  island,  total  surface  energy  of  the  system  and  their  changing  rates  associated  with  crack  growth  for 
(b)  l/l*r  =  1.1,  (c)  l/l*r  =  2/ \/3,  (d)  l/l*r  =  1.5;  (e)  phase  diagram  indicating  whether  a  preexisting  crack  will  grow  or  not,  and  if  so  whether  it  will  lead  to  full  delamination, 
for  given  crack  length  and  island  size. 


indicating  that  no  center  crack  can  grow,  regardless  of  its  length, 
when  the  island  size  is  small  enough.  For  the  case  of  /  >  l*r,  a  pre¬ 
existing  crack  of  length  a0  max  >  a0  >  a0  min  will  grow.  The  crack  will 
arrest  when 


Uel(a  =  a0)  -  l/e/ (a  =  Clan)  =  Usurf(a  =  aan)  -  Usurf(a  =  a0),  (26) 

where 


Usurf  =  r(a  -  a0),  (27) 

is  the  surface  energy  of  the  system  with  respect  to  the  initial  state. 
Eq.  (26)  can  be  solved  for  aarr  to  obtain 

aan  =  \  ^3/  -  2a0  +  J9P  +  12a0(  -12 ag  -  .  (28) 


By  comparing  aarr  with  l ,  we  can  determine  the  condition  for  com¬ 
plete  delamination.  Bearing  in  mind  the  condition  for  crack  growth, 
ao,min  <aoK  ao,max>  complete  delamination  aarr  >  l  requires 

^0,min  <  Oq  <  Qo,de/*  (29) 

where 


1 ,  /  /“  48  FEh  \ 

a°’de'-4  V  V  (1  -  V2)t^/2  J 

This  equation  breaks  down  for 


(30) 


r> 


3(1  -  v2)/2Tq 


16  Eh 


or  '  <  'a. del 


(31) 


indicating  that  full  delamination  is  impossible  for  sufficiently  small 
islands.  As  expected,  Eq.  (25),  which  inhibits  any  crack  growth, 
gives  a  stronger  condition  than  Eq.  (31),  which  allows  growth 
but  prevents  full  delamination.  Fig.  8a  shows  variations  of  the 
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Table  1 

Numerical  estimates  of  the  minimum  length  of  the  pre-existing  crack  required  for 
full  delamination  for  different  island  sizes  based  on  the  plane  strain  analysis. 


l/l*  „  , 

1  cr,del 

1 

1.34 

1.67 

2 

Island  size 

21  (|xm) 

29.95 

40 

50 

60 

^0,min  (p-ITl) 

3.74 

2.39 

1.81 

1.47 

Max.  stress 

crmax  (GPa) 

1.12 

1.76 

2.32 

2.85 

dimensionless  elastic  energy  release  rate  as  a  function  of  the 
dimensionless  crack  length.  As  shown  in  this  figure,  the  elas¬ 
tic  energy  release  rate  falls  entirely  below  the  rate  at  which  the 
surface  energy  of  the  system  increases  for  l  <  l*r.  However,  for 
l  <  l*r ,  these  two  curves  intersect  each  other  at  two  points  cor¬ 
responding  to  a0>min  and  a0,max-  The  total  elastic  energy,  the  total 
surface  energy,  and  their  changing  rate  as  a  function  of  crack 
length  are  shown  for  various  island  sizes  in  Figs.  7b-d.  In  each 
of  these  figures,  the  surface  energy  curve  has  been  plotted  for 
two  different  preexisting  crack  lengths  a0  =  a0m[n,  and  a0  =  0.5/.  As 
shown  in  Fig.  8b  for  l  =  1.1  Z*r  <  l*  d  r  neither  of  these  two  pre¬ 
existing  cracks  can  grow  all  the  way  to  the  edge,  as  the  elastic 
energy  in  the  patch  is  not  sufficient  to  drive  full  delamination, 
i.e.  Uei(a  =  do)  -  Uei(a  =  l)<  Usurj{a  =  Z)  -  Usurf(a  =  a0).This  is  consistent 
with  the  critical  condition  for  complete  delamination  given  by 
Eq.  (31).  In  the  critical  case  of  Z  =  (2/V3)l*r  =  l*r  deV  shown  on 
Fig.  8c,  a  preexisting  crack  with  a0  =  ao,min  =  0.25/  can  grow  all 
the  way  to  the  edge  leading  to  complete  delamination,  as  Eq. 
(26)  is  satisfied.  In  comparison,  a  preexisting  crack  with  a0  =  0.5Z 
is  arrested  when  LZe/(a  =  a0)-  Uei(a  =  l)<Usurj(a  =  l)-  Usurf(a  =  a0). 
Fig.  8d  shows  the  case  of  l  =  1 .5/Jr  >  Z*r  d  V  in  which  both 
of  the  preexisting  cracks  considered  here  grow  all  the  way 
to  the  edge,  leading  to  full  delamination,  as  they  satisfy 
Cl  —  CIq )  —  FZe/(fl  —  Z)  >  USUrf{ Q  ~  0  ~  ^surfi®  —  @0 )• 

A  summary  of  the  analysis  presented  above  can  be  cast  in  a 
phase-diagram  shown  in  Fig.  8e.  In  this  figure,  the  horizontal  axis  is 
a  measure  of  the  island  size  or  interface  toughness,  and  the  vertical 
axis  shows  the  length  of  the  preexisting  crack.  Three  regions  have 
been  indicated  on  this  figure,  which  are  separated  by  solid  curves. 
Given  the  dimensionless  length  of  the  preexisting  crack  a0/Z  and  the 
dimensionless  size  of  the  island  Z*r/Z,  the  phase  diagram  indicates 
whether  the  crack  can  grow  or  not,  and  if  it  can,  whether  it  will  be 
arrested.  The  dashed  curve  labeled  as  “2a0  min/Z”  provides  an  exam¬ 
ple  of  the  case  of  a0  =  2a0  m[n.  As  shown  in  the  figure,  this  curve  goes 
through  all  three  regions.  For  the  segment  of  this  curve  which  falls 
in  the  “Growth  &  Arrest”  region,  the  length  of  the  arrested  crack 
has  been  shown  by  the  dashed  curve  labeled  as  “(aflrr/Z)”,  which  can 
be  read  out  from  the  vertical  axis  for  a  given  Z*r/Z. 

For  lOOnm  thick  Si  islands  with  F  =  40GPa  and  v  =  0.22  [32], 
To  « 10  MPa  [21,22]  and  r  =  1  J/m2,  the  critical  island  sizes  2Z*r  and 
2 1*  dd  are  estimated  to  be  25.94  p,m  and  29.94  |jim,  respectively. 
The  minimum  length  of  the  required  pre-existing  cracks  for  full 
delamination,  and  the  corresponding  maximum  stresses  in  the  film 
prior  to  crack  growth,  are  given  for  different  island  sizes  in  Table  1 .  It 
should  be  noted  that  while  the  maximum  stresses  given  in  this  table 
exceed  the  apparent  flow  stress  observed  in  Si  thin  films  (1-1.7  GPa) 
[33],  our  assumption  of  overlapping  interfacial  sliding  in  the  island 
is  not  necessarily  violated.  This  is  because  the  latter  is,  in  fact,  an 
averaged  flow  stress  while  the  actual  stress  distribution  in  the  film 
is  not  uniform. 

4.2.  Axisymmetric  analysis 

In  this  section,  we  consider  a  central  penny-shaped  interfacial 
crack  of  radius  a,  as  shown  in  Fig.  9.  Invoking  the  same  assumption 
of  overlapping  interfacial  sliding  zones  as  in  the  foregoing  section, 


uniform  radial 
shear  stress 

Fig.  9.  Schematic  of  (a)  a  circular  Si  island  on  substrate  with  a  penny-shaped  center 
crack  and  (b)  a  free-body  diagram  of  the  island. 


Appendix  B  shows  that  radial  and  hoop  stresses  within  the  patch 
can  be  expressed  as 


CTrr  — 


TO 


6  hR2 


2(2  +  v)R3  -  3(1  +  v)R2a  -  a3{  1  -  v) 
(R-r)  [2(2  +  v)R2+(l-v)^-(R  +  r) 


a  >  r  >  0 
R> r > a ’ 


(32) 


Gee 


To 

6  hR2 


2(2  +  v)R3  -  3(1  +  v)R2a  -  a3(l  -  v)  a  >  r  >  0 

a3 

2(2  +  v)R3  -  2(1  +  2 v)rR2  -  (1  -  v)  —  (R2  +  r2)  R>  r  >  a 

r2 


(33) 


The  total  elastic  energy  of  the  patch  can  thus  be  calculated  as 


7l(  1  -  V)t3  -  ,  , 

Ue,  =  ((5  +  v)R6  -  8(2  +  v)R3a3  +  9(1  +  v)R2a 4  +  2(1  -  v)a6). 

36 EhR2 

The  elastic  energy  release  rate  is  then 


(34) 


G  = 


du ; 


el 


jr(  1  —  V)Tq  2 


da 


3  EhR2 


az{2R3{2  +  v)  -  3(1  +  v)R2a  -  (1  -  v)a3). 

(35) 


During  crack  growth,  the  surface  energy  associated  with  the  crack, 
Usmf=n{a2-al)r,  (36) 

would  grow  at  a  rate  of 

d^L  =  2  nar.  (37) 

da 

While  the  elastic  energy  release  rate  and  its  first  derivative  with 
respect  to  the  crack  length  vanish  at  a  =  0,  dUsurflda  vanishes  at 
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Fig.  10.  Axisymmetric  analysis  of  interfacial  delamination  of  a  circular  island  on 
substrate,  (a)  Critical  length  scales  as  a  function  of  the  Poisson’s  ratio  of  the  island; 
(b)  phase  diagram  in  terms  of  preexisting  crack  length  and  island  size. 


a  =  0  but  has  a  finite  positive  derivative  with  respect  to  a.  There¬ 
fore  G<dUsurflda  in  the  vicinity  of  a  =  0,  and  similar  to  the  plane 
strain  case,  cracking  at  the  center  is  impossible  without  a  preexist¬ 
ing  crack.  While  essential  features  of  this  problem  remain  similar 
to  the  plane  strain  analysis  presented  before,  the  algebraic  calcula¬ 
tions  become  rather  complicated  since  the  elastic  energy  is  a  sixth 
order  polynomial  with  respect  to  the  crack  length.  The  details  of 
the  analysis  are  skipped  here  and  the  emphasis  is  directed  at  the 
main  findings.  Similar  to  the  plane  strain  case,  we  can  identify  two 
critical  length  scales  R*r,  and  R*r  ...  If  the  island  radius  falls  below 
R*r ,  no  preexisting  circular  crack  can  grow,  regardless  of  its  radius, 
while  if  the  island  radius  falls  below  R*r  d  V  no  preexisting  crack  can 
lead  to  complete  delamination.  These  length  scales  as  a  function  of 
the  Poisson’s  ratio  of  the  island  are  shown  in  Fig.  10a.  For  lOOnm 
thick  Si  islands  with  E  =  40  GPa  and  v  =  0.22  [32],  ro  ~  10  MPa  [21,22] 
and  r  =  1  J/m2,  numerical  estimates  for  the  critical  island  sizes  are 
2  R*r  ^31.2  p>m  and  2  R*r  dd  %  37.6  p>m.  The  minimum  length  of  the 
center  cracks  which  leads  to  complete  delamination,  and  the  corre¬ 
sponding  maximum  stresses  for  different  island  sizes,  are  given  in 
Table  2.  A  summary  of  the  detailed  analysis  is  presented  in  Fig.  10b, 
in  the  form  of  a  phase  diagram  similar  to  the  one  presented  for  the 
plane  strain  case,  for  the  case  v  =  0.25. 


Table  2 

Numerical  estimates  of  the  minimum  length  of  the  pre-existing  crack  required  for 
full  delamination  for  circular  islands  of  different  sizes. 


R/Kr,del 

1 

1.06 

2.12 

3.19 

Island  size 

2 R  (|mm) 

37.64 

40 

60 

80 

^O.min  (l^ni) 

4.63 

4.20 

2.48 

1.80 

Max.  stress 

ermax  (GPa) 

1.11 

1.22 

2.06 

2.85 

5.  Conclusion 

In  summary,  we  have  proposed  a  method  to  deduce  the  critical 
size  for  interfacial  delamination  of  patterned  thin  film  electrodes 
on  substrate  in  the  presence  of  large  scale  interfacial  sliding.  We 
have  investigated  the  critical  conditions  of  interfacial  delamination 
with  plane  strain  and  axisymmetric  analysis  of  both  edge  and  center 
delamination  cracks  at  the  interface  between  thin  film  islands  and 
the  substrate.  Our  analysis  shows  that  the  overlapping  interfacial 
sliding  zones  due  to  the  finite  interfacial  sliding  strength  lead  to 
a  new  class  of  size  effects  in  interfacial  delamination.  Interfacial 
delamination  of  the  electrodes  becomes  impossible  when  the  island 
size  falls  below  a  critical  length  scale.  Applications  to  lithiation  of 
thin-film  Si  islands  give  results  in  excellent  agreement  with  the 
experimentally  observed  size  effects.  The  present  work  provides 
insights  for  optimized  structural  design  to  mitigate  the  mechanical 
degradation  in  high-capacity  electrodes. 
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Appendix  A. 


Consider  a  circular  membrane  of  radius  Ra ,  subject  to  a  uniform 
mismatch  strain  s *  and  uniform  radial  shear  stress  To.  Deformation 
of  the  patch,  within  plane  stress  approximation  and  in  the  pres¬ 
ence  of  interfacial  shear  traction  as  body-force  per  unit  volume  of 
the  patch  b  =  r0//i,  can  be  described  by  the  following  differential 
equation 


d_ 

dr 


i  _a 

r  dr 


(ru) 


+ 


b{\  -  v2) 


=  0, 


(Al) 


where  u  is  the  radial  displacement.  General  solution  to  the  differ¬ 
ential  Eq.  (Al )  can  be  written  in  the  form: 


u  =Ar  + 


B 

r 


1  Tq(1  —  V2)  2 
3  Eh 


(A2) 


where  A  and  B  are  constants  to  be  determined  from  boundary  con¬ 
ditions.  Radial  and  hoop  stresses  in  the  patch  are  related  to  the 
displacement  field  via 


Gnr 


E 

1  -  v2 


du  u 

—  +  V - (1  +  v)£* 

dr  r 


(A3) 


a00 


E 

1  -  V2 


du  u 

v -5-  +  -  -(1  +  v)e* 
dr  r 


(A4) 


Upon  satisfying  the  boundary  conditions,  u  =  0  at  r=0,  and 
<rrr  =  0  at  r=Ra ,  the  constants  A  and  B  are  obtained  as 
A  =  s*  +  ( 1  /3)r 0Ra(  1  -  v)(2  +  v)/E/i,  B  =  0,  and  subsequently,  the  stress 
components  are  found  as 

CT(t=  ro(2  +  y)(Ra_r);  (A5) 

CT00  =  ^i((2  +  v)Ra-(1+2v)r)).  (A6) 

Finally,  Eqs.  (10)  and  (11)  can  be  reproduced  by  setting  Ra  =  R-a. 


Appendix  B. 

Consider  a  membrane  of  radius  R,  subject  to  a  uniform  shear 
stress  to  within  the  region  R>r>a  and  zero  shear  stress  within 
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a>r> 0.  In  this  case,  the  body-force  in  Eq.  (Al)  should  be  replaced 
by 


b  = 


0  a  >  r  >  0 

^  R> r >  a 
h 


(Bl) 


This  equation  can  be  solved  within  R  >  r  >  0  to  give 

a  >  r  >  0 


i/  = 


A\r  +  — 
r 


„  B2  1  To(l  -  V2)  2  „ 

A2r  +  —  -  -  - -r2  R>r>a 

r  3  Eh 


(B2) 


where  At ,  A2 ,  £>i  and  B2  are  constants  to  be  determined  from  bound¬ 
ary  conditions  u  =  0  at  r  =  0,  and  <rrr  =  0  at  r=fta,  and  continuity 
conditions  of  the  displacement  as  well  as  the  radial  stress  at  r  =  a. 
This  yields 

^  =  e*  -  lTo{] L~  V)(2R3(2  +  v)  -  3aK2(l  +  v)  -  a3(l  -  v)), 

6  Eh 


A2  =  A\  + 


r0a(l  -  v2) 


B\  =  0  and  B2  = 


r0a3(l  -  v2) 


2Eh  6 Eh 

Upon  substitution  of  these  constants  in  Eq.  (B2),  the  radial  and  hoop 
stresses  can  be  obtained  from  Eqs.  (A3)  and  (A4),  as  given  in  Eqs. 
(32)  and  (33). 
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